Semiclassical Approach to Survival Probability at Quantum Phase Transitions 



Wen-ge Wang*, Pinquan Qin, Lewei He, and Ping Wang 
Department of Modern Physics, University of Science and Technology of China, Hefei, 230026, China 

(Dated: February 2, 2010) 

We study the decay of survival probability at quantum phase transitions (QPT) with infinitely- 
degenerate ground levels at critical points. For relatively long times, the semiclassical theory predicts 
power law decay of the survival probability in systems with d — 1 and exponential decay in systems 
with sufficiently large d, where d is the degrees of freedom of the classical counterpart of the system. 
The predictions are checked numerically in four models. 
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I. INTRODUCTION 

A quantum phase transition (QPT) is characterized by 
non-analyticity of the ground level of the system at the 
critical point in the large size limit. At a QPT, certain 
fundamental properties of the ground state (GS) change 
drastically under small variation of a controlling param- 
eter, e.g., strength of a magnetic field. Most of the works 
in QPT have focused on properties of equilibrium states 
(including GS at zero temperature) [1]. While, the non- 
analyticity influences in fact both equilibrium and non- 
equilibrium properties. Indeed, when the time scale of in- 
terest is smaller than the relaxation time, which diverges 
at the critical point, usually the system is not in an equi- 
librium state and unitary dynamics should be considered 
(FigO}. Due to significant progress in cold atom experi- 
ments, time dependent simulation of models undergoing 
QPT is becoming realizable 0, H[, hence, investigation 
in the unitary dynamics at QPT is of interest both the- 
oretically and experimentally. For example, resorting to 
theoretical technique such as a quantum version of the 
Kibble-Zurek theory 0, 0], it has been shown that slow 
change of the controlling parameter passing the critical 
point may induce some intriguing effects Q. 

In this paper, we study a different dynamics at QPT, 
which is induced by a sudden small change in the con- 
trolling parameter, A — > A', in the vicinity of a critical 
point A c . A measure of the effect of this dynamics is the 
survival probability (SP) of an initial state prepared in 
the GS |0a) of H(X), 

M(t) = |(0 A | e -^ A '>*/ fi |0 A )| 2 . (1) 

The SP, sometimes called autocorrelation function, is 
a quantity accessible experimentally Q. Recent it was 
found that relatively significant and fast decay of the SP 
may indicate the position of QPT [8rJIo| , which has been 
demonstrated experimentally [TTI |. Short time decay of 
the SP has been studied in these works. Of further in- 
terest, while still unknown, is the law for relatively- long- 
time decay of the SP at QPT and whether it may be use- 
ful in revealing characteristic properties of QPT [l2|, [HI ■ 
To find an answer to the above question, here we focus 
on those QPT, at the critical points of which the ground 
levels have infinite degeneracy in the large size limit. This 
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FIG. 1: A schematic plot, where r is the time scale of interest 
and T r is the relaxation time, A is a controlling parameter with 
critical value A c of a QPT. Below the solid curves, r < r r , the 
system is usually not in an equilibrium state and its unitary 
dynamics should be considered. 

is a type of QPT met in many cases (see models discussed 
below and those in Ref. |l|). At such a QPT, the non- 
analyticity may be a consequence of avoided crossings of 
infinite levels, not a few levels. 

We find that the semiclassical theory may be used in 
the study of the SP decay when A' is sufficiently close 
to A c . The theory predicts a power law decay of the 
SP in some systems and an exponential decay in some 
other systems. Numerical results obtained in four models 
confirm these predictions. 



II. SEMICLASSICAL APPROACH 

We first discuss a condition for the applicability of the 
semiclassical theory in the study of the SP of GS. We 
use notations: e = A' — A, 6 = A' — A C ,AA = A — A c 
(FigHJ, and r\ = e/AA. We use \et\) with a = 0, 1, . . . 
to denote eigenstates of H(X) with eigenenergies E a (X) 
in increasing energy order. When the ground level of 
H(X C ) is infinitely degenerate and those of H(X') are non- 
degenerate (or have finite degeneracy), infinitely many 
low-lying levels of H(X') must join its ground level in the 
limit A' — > A c , i.e., 

lim E a (X') = Eq(X c ), for many a. (2) 

This has two consequences: (i) H(X') of A' sufficiently 
close to A c must have a high density of states near its 
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ground level, (ii) For a fixed A near A c , when A' is suffi- 
ciently close to A c , H(X') may have many levels below E, 
where E = {Q\\H{\')\Q\), i.e., the initial state |0x) may 
have a relatively high mean energy in the system H(X'). 
This is in agreement with a property revealed in recent 
study of the fidelity of GS near critical points, which has 
close relationship to the SP, namely, for a fixed small 
e, the overlap |(CK|0a')I decreases significantly when A' 
approaches A c [14j . 

Moreover, suppose the system has a classical counter- 
part in the low energy region. Here, a classical counter- 
part means a classical system, the quantization of which 
gives a system mathematically equivalent to the origi- 
nal quantum system; its components are not required to 
be directly related to components of the original system. 
The property ([2]) implies that in the process A' — > A c 
longer and longer trajectories in the classical system may 
be of relevance. For a fixed initial state \0\), one may as- 
sume that the initial value of the Lagrangian L does not 
change notably in this process. Then, trajectories of rele- 
vance may have large action S — J Ldt' for A' sufficiently 
close to A c . 

The above discussed properties for A' sufficiently close 
to A c , namely, high density of states, relative highness of 
E, and large action of some relevant classical trajectories, 
imply that a semiclassical approach may be valid. To be 
specific, for any given A near A c , it is reasonable to expect 
that the semiclassical theory may be used in the study of 
the SP when A' is sufficiently close to A c . 

According to the semiclassical theory, qualitative dif- 
ference in classical trajectories may have quantum man- 
ifestation. Specifically, in the case of d = 1 where d is 
the degree(s) of freedom of the classical counterpart in 
the configuration space, the classical motion may show 
periodicity within a time scale of interest; on the other 
hand, for a large d, even in a regular system, classical 
trajectories may show no signature of periodicity within 
times of practical interest. This difference suggests that 
the SP decay in the former case may be slower than in 
the latter case, which we discuss below. 

We consider small e, such that H(X') = H(X) + eV, 

with V ~ ^§£^. The SP of the GS of H(X) is a special 
case of the so-called quantum Loschmidt echo or (Peres) 
fidelity [7J], M L (i) = |m(t)| 2 , where 

m(t) = {^ \exp(iH(X')t/?i)exp(-iH(X)t/h)\^o}- (3) 

In studying the SP, one may employ a semiclassical ap- 
proach that has been found successful in the study of 
Loschmidt echo [l6r 23]. For an initial Gaussian wave 
packet, narrow in the coordinate space with width a and 
centered at (?o,po) in the phase space, using the semi- 
classical Van Vleck-Gutzwiller propagator, it has been 
shown that [H [TJ] 
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for small e, which works in both regular and chaotic cases 
18, 23]. Here, AS is the action difference between two 



nearby trajectories in the two systems starting at (po, r ) 
and approximately can be evaluated along one trajectory, 
AS ~ e f*dt'V[r(t'),p(t')} pjj. The quantity w is h/a 
for sufficiently small a and depends on both a and the 
local instability of the classical trajectory when a is not 
very small [201 ] . 

We first discuss the SP in the case of d = 1 with a reg- 
ular dynamics. We assume that the GS can be (approx- 
imately) written as a Gaussian wave packet in certain 
coordinate of the classical counterpart. This is possi- 
ble, e.g., in the models discussed below. In this case, as 
shown in Ref. [HI, for t > T, due to the periodicity of the 
classical motion, the main contribution of AS to the SP 

is given by its average part eUt, where U — ^ J Q T V(t)dt 
and T is the period of the classical motion in H(X). Upto 
the first order expansion of U in p n , Eq. @ predicts a 
Gaussian decay of the SP [2l|, [23|. For relatively long 
times, higher order terms of U induces power law decay 
of the SP [23, 24]. For example, to the second order term, 



-i/2 e -rt 2 /(i-K 2 * 2 ) 



(5) 



where cq ~ 1, T 
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It is seen that Mi has 



indicating evaluation at po 2 
a Gaussian decay e _rt for initial times and has a l/£t 
decay for long times. 

Next, we consider the case of a regular classical coun- 
terpart with large d. In this case, the underlying classical 
motion is typically quasi-periodic with many different fre- 
quencies, as a result, T is usually much longer than time 
scales of practical interest. For times t -C T, classical 
trajectories may look random in the torus, due to the 
difference in the frequencies. To calculate the SP in this 
case, one may write it in terms of the distribution P(AS) 
of AS (with the Gaussian weight taken into account), 

M sc (t) ~ |J dASe iAS / h P{AS)\ 2 . When the trajectories 
can be effectively regarded as random walks for times 
f < T due to the many frequencies, P(AS) is close to 
a Gaussian distribution, independent of the initial state. 
In this case, the SP can be calculated in the same way as 
in a chaotic system (l7| , which has an exponential decay 
determined by the variance of AS*, 



M 2 (t) 



-K s e 2 t/h 2 



where 




Vdt 




(6) 



(7) 



with / Vdt = f*dt'V[r(t'),p(t')] [H. 

To summarize, for small e and sufficiently small 5, and 
for relatively long times, the SP may have a power law 
decay when d= 1, and has the exponential decay M2(t) 
when d is sufficiently large. We remark that, for A far 
from A c , the SP is always close to 1 for small e. 
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FIG. 2: (Color online) Decay of the SP (dashed curves) in 
the normal phase of Dicke model. Parameters: u> = ujq = 1, 
e = 1CT 5 , and 5 = -l(T m with m = 6, 7, 8, 9, 10, 11 from top 
to bottom. The solid curve is a fitting curve of the form in 
Eq. (0, having an initial Gaussian decay e~ Ft followed by a 
l/£t decay. The l/t decay becomes clear with increasing m, 
i.e., with A' approaching A c . Upper right inset: (In M)/et 2 
versus r\ for different pairs of (e, i) with short t, in agreement 
with the prediction P ~ \rje\. Lower left inset: InM versus 
ln(e 1/2 i) for e G (10" 6 ,10 -5 ) and lnt € (8.6,9.5) in the l/t 
decay region. 8 — — 10~ 10 , thus, \r/\ ~ 1. The results are in 
agreement with the prediction £ ~ jryel 1 ' 2 . 



III. NUMERICAL SIMULATIONS 

The first model we study is the single-mode Dicke 
model 26], describing the interaction between a single 
bosonic mode and a collection of N two-level atoms. In 
terms of collective operators J for the N atoms, the Dicke 
Hamiltonian is written as (hereafter we take fi = 1) [27| . 

H = lo q J z + Loa f a + (X/VN)(a 1 + o)(J+ + J_). (8) 

In the limit N — > oo, the system undergoes a QPT at 
A c = ^^/ujujq, with a normal phase for A < A c and a 
super-radiant phase for A > A c . The Hamiltonian can be 
diagonalized in this limit, 

H ( x ) = e ^ c lx c kx + g, (9) 

fe=l,2 

where ct\ and c k x are bosonic creation and annihilation 
operators, e k x are single quasi-particle energies, and g is 
a c-number function [27] . To be specific, in the normal 
phase, 

elx = \ + + (-l)V^o-^ 2 ) 2 + 16A 2 ^ | . 

(10) 

It is seen that e\\ c = 0, hence, the ground level of H(X C ) 
is infinitely degenerate. Since e2x c = \/u 2 + oJq is finite, 
at the QPT one may consider the effective Hamiltonian 
Hes(X) = eix c \\Cix with d = 1. Direct calculation shows 



e u =s AIAAI 1 / 2 , with A = 2( ;" o)3/ 2 4 , and 

y = - 2^: ( ct ^ ciA + 2(ct ^ )2 + 2c? ~ |AArl/2 - 

(ii) 

The semiclassical result Eq. ([5]) predicts that the SP 
has a Gaussian decay followed by a power law decay, with 
scaling properties 

r ~^p = M> £~M 1/2 . (i 2 ) 

These predictions have been confirmed in our numeri- 
cal simulations (Fig. [2]). Numerically, the SP was calcu- 
lated by making use of relations between (c|, A , c^) and 
(ct v , Cfex' ) i which can be directly derived from formu- 
las given in Ref. (27l |. Our numerical results support the 
prediction that the semiclassical theory may work for suf- 
ficiently small \S\. Similar results have also been found 
in the super-radiant phase. 

The second model we have studied is the LMG model 
[H], with the Hamiltonian H = -jASl + 7^) - XS Z , 
which has a critical point at A c = 1[22]. The model has 
a classical counterpart with d = 1. Direct computation 
shows a l/t decay of the SP for relatively long times in 
the neighborhood of A c [30]. 

As a third model, we study a 1-dimensional Ising chain 
in a transverse field, 

JV 

H(X) = -J2a!** +1 +Xaf. (13) 

i=i 

The Hamiltonian can be diagonalized by using Jordan- 
Wigner and Bogoliubov transformations, giving H(X) = 

J2k ek (^k^ k ~ 1/2) [!]■ Here, b\ and are creation 
and annihilation operators for fermions and are single 
quasi-particle energies, 

e k = 2Vl + A 2 -2Acos(fca) (14) 

with lattice spacing a, where k = with m = 

—M, -M + l,...,M and N = 2M + l.°Note that (ka) 
in Eq. (fT4|) is in fact independent of the lattice spacing 
a, with ka = 2-nm/N . 

To understand the degeneracy property of the ground 
level in the large N limit, let us consider those m satis- 
fying |m| < N@ for large TV, where /3 € (0, 1) is an arbi- 
trary number independent of N . In the limit N — > oo, 
one has ka — > for these m. As a result, Eq. (fT4")) gives 
efc = 2|AA| with A c = 1, in particular, at the critical 
point A = A c , efe = for these modes m. The number 
of these modes m is infinitely large in the limit f->oo, 
hence, the ground level is infinitely degenerate. 

In a sufficiently low energy region and for A ~ A c , 
a classical counterpart of the system can be introduced 
as follows. For A = A c , e k — 4:ir\m\/N for sufficiently 
large N and small \m\. In the low energy region, due 
to this linear dependence of e k on m, using the method 
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where 



FIG. 3: (Color online) Decay of the SP (circles) in a 1- 
dimensional Ising chain in a transverse held, with N = 2 X 10 s , 
e = 8x 10 -6 , and 8 = — 4x 10 -6 . It has the expected exponen- 
tial decay (solid straight line). Lower left inset: Dependence 
of ln(— InM) on ln|e| for a fixed time t. The straight line 
has a slope 2, as predicted in Eq. ([6]). Upper right inset: The 
SP increases slowly with increasing |(5| for fixed e and t, in 
agreement with the prediction for K s given in the text. 



F k = 1 - sirT 



= arctan 



A') sin (e fe t), 
— sin(fca) 



(18) 
(19) 



cos(fca) — A ' 

and e k are evaluated at A' . Our numerical computations 
confirm not only the prediction of an exponential decay 
of the SP at the criticality, but also some details in the 
exponent of M 2 (t) discussed above (see FigJ3]), namely, 
the e 2 dependence and the properties of K s . 

As a fourth model, we have studied the XY model [l|, 
with the Hamiltonian 



1 + 7 



1 



'i+l 



> V V 



2 Ui 



(20) 



which has critical points A c = ±1. As in the Ising chain, 
in the low energy region around A c , the XY model has a 
classical counterpart with large d. The SP in this model 
can be calculated in a way similar to that in the Ising 
chain discussed above, and our numerical simulations also 
confirmed the semiclassically predicted exponential decay 
of the SP. 



of bosonization (see Ref. [l|), one can express fermionic 
states ¥ ki . . . b' k | vacuum) in terms of (many) bosonic 
modes. Each bosonic mode has a classical counterpart 
with one degree of freedom, hence, H(X C ) has a classical 
counterpart in the low energy region with a large value of 
d (d — > oo in the large N limit). This implies that H{\) 
with A ~ A c also has a classical counterpart with large d, 
as a result, typically the SP should have an exponential 
decay M 2 (t) in Eq. ©. 

Direct derivation shows that the perturbation in this 
model is 

A — cosfca., + , 1, sinfca.,, , jtlt 
V = — —— (blb k --) + ——i(b k b_ k -blbl k ). (15) 
efe/4 i ekjl 

Further analysis shows that V has no singularity at the 
critical point, e.g., 



sinfca J sinfca/|AA|, 



for \ka\ < |AA| 

sinfca/\/l — cosfca, for \ka\ > |AA| 



.(16) 



Therefore, K s in Eq. has no singularity in the vicinity 
of A c . For large and fixed N and for |AA| 1/N, since 
the coupling strength of V in the eigenbasis of H (A) in- 
creases with decreasing |AA|, K s should increase slowly 
with decreasing |AA|. 

Numerical computation of the SP can be done by using 
the following expression given in Ref. @ , 



m (*) = n F k 



(17) 



IV. CONCLUSIONS AND DISCUSSIONS 

We have shown that the semiclassical theory may be 
used in the study of the decay of SP (survival probabil- 
ity) of GS (ground states) in the vicinity of those QPT 
with infinitely degenerate ground levels at the critical 
points. Two qualitatively different decaying behaviors of 
the SP have been found for relatively long times: power 
law decay in systems with d = 1 and exponential decay in 
systems with sufficiently large d, where d is the degrees 
of freedom of the classical counterpart of the quantum 
system. 

The above results suggest that the SP decay may 
be useful in the classification of QPT, an important 
topic far from being completely solved, in particular, in 
the non-equilibrium regime. Here, we have found two 
classes: one class with power law decay and another class 
with exponential decay. It needs further investigation 
whether other types of SP decay may appear at QPT, 
e.g., relatively- long-time Gaussian decay or a decay be- 
tween power-law and exponential. 
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